Quantum transport anomalies in DNA containing mispairs 
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Abstract 

The effect of mispair on charge transport in a DNA of sequence (GC)(TA)7v(GC)3 connected 
to platinum electrodes is studied using the tight-binding model. With parameters derived from 
ah initio density functional result, we calculate the current versus bias voltage for DNA with and 
without mispair and for different numbers of (TA) basepairs between the single and triple (GC) 
basepairs. The current decays exponentially with N under low bias but reaches a minimum under 
high bias when a multichannel transport mechanism is established. A (GA) mispair substituting 
a (TA) basepair near the middle of the (TA) tv sequence usually enhances the current by one order 
due to its low ionization energy but may decrease the current significantly when an established 
multichannel mechanism is broken. 
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I. INTRODUCTION 



Longitudinal charge transport along DNA has been the subject of extensive study in the 
last decade, [l-st Charge transport occurs in the oxidative and reductive DNA damage 
or repair processes and can happen in the long distance range. Study of trans- 

port properties may lead to a better understanding of the fundamental driving processes 
in biological evolution. Furthermore, the charge transport process might have been used 
naturally for basepair mismatch detection during the DNA repairing process. It is already 
known that, due to chemical reaction and radiative ionization, mispairs or gene mutations 
happen quite often in the cell. Fortunately, almost all of the mispairs can be detected and 
repaired during the replication process to keep the material genetically stable. However, 
some of the mutations may escape from the detecting and repairing processes and result in 
various genetic diseases including cancer. A recent study indicates a negative correlation 
between the cancer risk and sensitivity of charge transport property of the gene to a muta- 
tion. [6] Understanding how mispairs modify the electric properties of DNA then becomes 
very important 7-^ and, together with the usage of other properties, jiol . ll| may improve 



mutation detecting techniques. 
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In addition, thanks to its perfect self-assembling and 
self-recognition properties found in nature, DNA is also expected to be a potentially func- 
tional material for molecular devices. In this case mispairs may be used to obtain unique 
functions of the devices. 

The charge transport through a DNA sequence can be measured by chemical or physical 
methods. lH3|] In one of the typical chemical experiments, Giese et al. used a DNA of 
sequence (GC)(TA)Af(GC)3. |l6| They measured the charge transfer rate from the (GC) 
basepair to the (GC)3 triple basepair for different number of (TA) basepairs, and found 
a crossover from a rapid decay of the charge transfer rate vs N to an almost zero decay 
around = 3. As an alternative to other explanations, 17H22| we have proposed this as a 



crossover from one dominant channel transport to a multichanne 
of physical experiments is the one performed by Porath et al. 



transport. 23|] An example 



24| where a DNA sequence 



(GC)m is located between two platinum electrodes and the current versus voltage is directly 

hi It is 



measured. This result has also been simulated by simple tight-binding models. |25 
known that G.A mispair in various conformations ll| is the most stable mispair and often 
present in the DNA. js] The magnetic properties of DNA was studied earlier and found to 




FIG. 1: (a) Schematic illustration of the equilibrium energy band across the system with a DNA 
of sequence (GC)(TA)Ar(GC)3 connect to two platinum electrodes, (b) The energy dependence of 
em (solid curve) and tm (dashed) for electrons in the platinum electrode X, with X = L or R foi 
the left or right electrode. 



be significantly influenced by the presence of the G.A mispair. 10| In this paper, we will 
study the effect of mispairs, such as G(anti)-A(anti) indicated in the following as (GaAa), 
and G(anti)-A(syn) indicated as (GaAs), 28[ on charge transport when a Watson-Crick 
(TA) basepair is replaced by a mispair in the DNA sequence (GC)(TA)jv(GC)3 connected 
to platinum electrodes. 



II. METHOD 



We consider a p-type semiconductor DNA duplex chain of basepairs connected to a cir- 
cuit via two platinum electrodes suitable for experimental realization 24|. Each platinum 
electrode is modeled as a semi-infinite one-dimensional (ID) electrode 27| connected to the 
G base at one end of the first strand as illustrated in[l](a). The tight-binding Hamiltonian 
of the system reads 



H — 2 ^ ^ l^nCli^Cn — i^n,n+l(c]^C„+l + c|^+lCn)] 



N N-l 



n=l n=l 

N 



- 2^ A„(4ci„ + 4c„). (1) 



n=l 



Here cj^ (dj^) is the creation operator of holes in the first (second) strand on site n of 
the DNA chain (for 1 < n < N), the left electrodes (n < 0), and the right electrodes 
(n > + 1). The on-site energy of site n in the first (second) strand is denoted by (wn), 
which is equal to the highest occupied molecular orbit (HOMO) energy of the base on this 
site in the DNA chain and the center of conduction band in the electrodes. The coupling 
parameter of the first (second) strand t„,„+i (/in,n+i) is equal to the intra-strand coupling 
parameter between neighboring sites n and n + 1 of the DNA for 1 < n < A^ — 1 , one- fourth 
of the conduction band-width in the electrodes tm for n < —1 and n > N + 1, and the 
coupling strength between the electrodes and the DNA strands for n = and n = N. The 
inter-strand coupling between sites in the same basepair is described by A„. The factor 2 
multiplied to each sum in E g. Q arises from the spin degeneracy. 

In transport experiments, 2J] a high bias voltage can be applied to drive the system far 
from equilibrium and holes in wide energy range may contribute to the current. Since the 
carriers usually come from various energy bands and the profile of band distribution is energy 
dependent, the effective parameters Em and tm, which are averages over the profiles, are then 
energy dependent. We assume that the parameters for the ID tight-binding model have a 
similar dependence on energy as in bulk platinum j3] and the dependence is extracted from 
its 3D band structure. Near the Fermi energy, there are six bands located approximately at 
—5.8, —4.7, —3.7, —2.2, —0.2, and 2.0eV above the Fermi energy with band width 1.9, 1.3, 
1.5, 3.1, 1.4, and 6.0 eV respectively. Using Lorentzian broadening, we can mimic the bulk 
DOS and extract the parameters Em and tm as shown in[T]^b). The parameters are then scaled 
to match the known values at the Fermi energy as was done in Ref. 27|]. For electrons at 
the Fermi energy, the on-site energy is = —0.33 eV with a coupling parameter = 0.55 
eV. As estimated from the experimental data 25|, l27|] the equilibrium Fermi energy is 1.73 
eV higher than the HOMO on-site energy of the G base when the (G-C) basepair makes 
contact with the platinum electrodes. Here we assume that the first DNA strand is coupled 
to the electrodes with a contact parameter of tdm = 0.1 eV while the second strand does not 
contact the electrodes directly. Note that our main result is not sensitive to the choice of 
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FIG. 2: Illustration of a three-basepair DNA used in the ADF program to obtain the tight-binding 
parameters. In the first (second) strand, the on-site energy of a HOMO orbital is e„ and the 
intrastrand coupling parameter between neighboring sites is tn,n+i (^n,n+i) with n the base index. 
The interstrand coupling parameter is denoted by A^. 



the electrodes and the contact parameters. 

The tight-binding parameters of DNA are estimated based on the HOMO energies of 
isolated nucleobases and the charge transfer integral between the HOMO orbitals calculated 
by the ab initio density functional method integrated in the ADF (Amsterdam Density 



Functional) program. 10 
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30| The on-site energies for bases G, C, T, and A are —9.40, 



— 10.27, —10.46, and —9.79, respectively. The hopping coupling parameters are listed in 
Table E 

The current / when a voltage bias V is applied over the two platinum electrodes is then 



evaluated by the transfer matrix method 



I2, 23, lai 



32| . For an open system, the secular 



equation is expressed as a group of equations of the form 

with \1'+ (^^) the wave function of the first (second) strand on site n. The wave functions 
of the sites n + 1 and n are related to those of the sites n and n — 1 by a transfer matrix M, 



n+1 



M 



n 



(2) 



TABLE I: The values of the intra- and interstrand tight-binding parameters (in eV) for different 
DNA sequences, where the middle pair (Xi-X2) can be G-C and T-A basepairs or G(anti)-A(anti) 



m 


5'-G-Xi-G-3' 
3'-C-X2-C-5' 


(Xi-Xs) 


*12 ^23 ^^2 ^12 ^23 


GC 
T-A 


0.133 0.133 0.028 0.14 0.14 
0.164 0.400 0.070 0.154 0.099 


5'-T-Xi-T-3' 
3'-A-X2-A-5' 


(Xi-Xs) 


*12 *23 -^2 ^12 ^23 


T-A 

Ga-Aa 
Ga-As 


0.330 0.330 0.070 0.011 0.011 
0.400 0.164 0.029 0.055 0.002 
0.250 0.167 0.057 0.207 0.027 



with 



M 



-Xrt 



-An {E-en) 



tn,n-\-l 










^ri,n+ 1 
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The transmission is then calculated by assuming the plane waves propagating in the elec- 
trodes for the holes = Ae*''^"'' + Be''^^'"' for n < and = Ce*'=^"° for n > + 1 
in the left and right electrodes, respectively. Expressing the output wave amplitude C in 
terms of the input wave amplitude A and the transmission, 

, , |Cp sin(A;/ja) 
^ ^ " \A\^sm{kLa)' 

The distance between two neighboring bases along any DNA strand is a = 3.4 A. The net 
current primarily comes from the hole transmission between the electrodes' Fermi energies 
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FIG. 3: (a) Current / versus bias voltage V of the DNA sequence (GC)(TA)Ar(GC)3. iV = 1 - 11 
for curves from the top. (b) Current I versus (TA) basepair number N at fixed bias voltage V . 
The value of V is indicated beside each curve. 

n 

and is calculated as [33i] 



Here the Fermi function is f^{E) = 1/ exp[(£' — Ep)/kBT] with X = L or R and the room 
temperature T = 300 K. When a bias voltage V is applied between the two electrodes, the 
left (right) Fermi energy is assumed Ep = V/2 {Ep = —V/2). 

III. RESULTS AND DISCUSSIONS 

The results for the current (/) versus voltage (V) for DNA of sequence (GC)(TA)Ar(GC)3 
with = 1,2, 11 are shown in[3t^a). Each curve except for = 1 has steps at voltages 
1.45, 1.64, 1.85, and 2.02 eV indicating that the transport channels are mainly formed at 
four threshold voltages related to the four kinds of bases. In the = 1 curve, the second 
step is split into two ai V = 1.62 and 1.7 V and there is no step at 2.02 V. In a log scale, 
the two lower-energy steps have almost the same height for all the curves while the height of 
the two higher-energy steps increases with the number of (TA) base-pairs. At a bias voltage 
lower than the first threshold, the DNA works as an energy barrier for electron transport 
since the Fermi energy of both electrodes is located between the HOMO and LUMO energy. 
At a bias near the first and the second threshold (1.45 < 1.8 eV), HOMO channels of the 
(GC) basepairs become available for charge transport but the (TA) basepairs behave as 
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energy barriers for transport. At a bias higher than the third threshold {V > 1.85 eV), 
transport channels of (TA) basepair also participate in the charge transport across the DNA 
molecule. Along the curve = 1, we can hardly see the third and the fourth steps. The 
addition of (TA) basepairs can establish a network of bases due to the interstrand coupling 
and introduces additional transport channels, as was reported in Ref.|23|. Consequently, 
we observe the height increase of the third and fourth steps in the log scale. The current 
enhancement due to additional channels may compete with the exponential current decay 
with the length of the molecule. This results in a current minimum for A^ > 5, as clearly 
shown in Wljo) where the current is plotted versus A^ at various V. At a bias less than 
1.8V, the current decays exponentially with A^ with almost the same exponent. For a bias 
higher than 1.8V, the exponent decreases with A^ until it is almost zero for N > 5 under 
the bias V = 2.1V or higher. This crossover from a rapid to almost zero decay of the charge 
transfer versus the (TA) basepair numbers has been observed in Ref. jl6| with the chemical 
method and can also be observed in physical experiments as described in Ref. {2^. Different 
from our previous simplified model [23|] where uniform parameters and virtual electrodes are 
assumed, here we employ a more realistic model with the tight-binding parameters of DNA 
and electrode extracted from ab initio calculations. In addition, a variable bias voltage is 
applied between the two electrodes to obtain the I-V curve. Note that the role of diagonal 
interstrand hopping is relevant to the electron transport in DNA and its inclusion in the 



calculation might shift the I-V curve but not the conclusion. [SJ, |35l] 

To estimate the effect of mispairs (GaAa) and (GaAs) on the charge transport, we replace 
one of the (TA) basepairs near the middle of the (TA)7v sequence by a mispair and calculate 
the corresponding I — V curve. In HJ^a) and (b), the result for the [(A^ + 2)/2]th (TA) 
basepair replaced by a (GaAa) mispair is shown and in (c) and (d) the [(A^ + l)/2]th basepair 
is replaced by a mispair. Here [R] means extracting the integer part of a real number R. 
Note that the curves of odd A^ are the same inlH^a) and (c), corresponding to equal number 
of (TA) basepairs to the left and right of the (GaAa) mispair. For even A^, there are one 
more (TA) to the left of the mispair in (a) and one less in (b). With the replacement of 
the mispair, the current is greatly enhanced after the first threshold voltage, indicating that 
the G base in the mispair works as a tunneling bridge. However, for A^ > 5, the current 
decreases with the mispair at a bias higher than the fourth threshold voltage (2 eV) as shown 
in 111(a) and (c). This happens because the mispair destroys the resonant transport network 




FIG. 4: (a) The I - V curves of DNA sequence (GC)(TA)^/2(GaAa)(TA)(;v-i)/2 (GC)3 for = 
3, 4, 11. Here N/2 and (A^ — 1)/2 take only the interger part of the value, (b) Current I versus the 
number A^ at several selected bias voltage V indicated by the values beside each curve. The I — V 
curves and I — N curves for DNA sequence (GC)(TA)(jv-i)/2(GaAa)(TA)^/2(GC)3 are plotted in 
(c) and (d). 

formed by the periodic (TA) basepairs series, i.e. the mispair works as an impurity. A weak 
dependent current appears only for higher A^ and at a lower current value, as shown in 
m^b) and (d) when both the (TA) sequences besides the mispair form resonant transport 
network. 

In the presentee of a mispair, some steps shift and extra steps appear along the I — V 
curves. For A^ = 3 or a DNA of sequence (GC)(TA)(GaAa)(TA)(GC)3, the simplest case 
with a mispair, three steps at V = 1.41, 1.55, and 1.68 V appear in the I — V curve and 
the current is enhanced by more than one order at high voltage with the substitution by 
a mispair. For A^ > 5, however, the first three steps shift to V = 1.45, 1.53, and 1.64 V 
with the last one decaying over A^. Compared to the case without a mispair, one extra step 
appears at = 1.53 V. In the bias range 1.45 V < V < 1.53 V, the curves for A^ > 4 
are almost equally separated in the vertical log scale, indicating an exponentially decaying 
current with A^ as also shown inlD^b) and (d). Furthermore, the mispair location also affects 
significantly the current. The A^ = 4 curve has steps at the same position as for A^ = 3 
in m^a) but as the A^ > 5 curves in Hl^b). In the range 1.58 V < ^ < 2.2V, the even A^ 
curves are near the curve of A^ + 1 in (a) but near the curve A^ — 1 in (b). The change is 



also represented by shifted steps along the V = 1.6 V curves when comparing |l](b) to (d). 
This observation suggests that the current in this bias range is mainly determined by the 
length of the (TA)„ sequence between the left single (GC) basepair and the (GaAa) mispair. 
Under stronger bias, the current at first decays exponentially with and then fluctuates 
near a value slightly below 1 nA. This long-DNA current limit is about five times smaller 
than that observed in the system without a mispair. 

When the (GaAa) mispair is replaced by a (GaAs) mispair, the I — V curves show fewer 
steps as illustrated in Fig.O For A^ = 3 there are three steps aX V = 1.45, 1.58, and 1.68V 
while for A^ > 5 there are only two steps at ^ = 1.45 and 1.61V below bias 2V. Similar to 
the case of (GaAa) mispair shown in HI the A^ = 4 curve also has steps at the position of 
the A^ = 3 curve in (a) but at the positions of the A^ > 5 curves in (b). In addition, the 
current is mainly limited by the number of (TA) basepairs between the left (GC) basepair 
and the (GaAs) mispair in the range 1.6V< V < 2.2V. For V > 2.2V, no step in curves 
of A^ < 7 suggests again that one transport channel dominates in the short (TA)„ DNA 
sequence. For large A^, a series of steps appear in the I — V curves and the current do not 
decrease monotonically with A^, indicating the enhancement of current due to the increase of 
transport channel number with A^ can compensate the current decay with the length of each 
channel. The current decays exponentially at < 1.6 V but decays with steps in the range 
1.6 < < 2.2 V. At V > 2.2 V the current decays in short DNA and then fluctuates when 
multichannel tunneling mechanism dominates in the long (TA)„ sequence. Overall a (GaAs) 
mispair substitution changes less the current than that of a (GaAa) mispair especially in 
high bias since the intrastrand coupling parameter of a (GaAs) mispair is closer to that of 
a (TA) basepair. 

IV. SUMMARY 

In summary, we have studied charge transport through DNA connected to two platinum 
electrodes and contains mispairs within a realistic tight-binding scheme. The energy de- 
pendent tight-binding parameters for the electrodes are obtained by fitting the density of 
states near the Fermi energy of the material. The parameters of DNA are derived from the 
ab initio density functional calculation of the coupling between HOMO states in neighbor 
bases. When a (TA) basepair in the (GC)(TA)7v(GC)3 sequence is replaced by (GaAa) or 




FIG. 5: The same as in|3]is plotted for basepair mismatch (GaAs). 

(GaAs) mispairs, the current is usually enhanced due to the lower ionization energy of the 
mispairs. In DNA with a long (TA)7v sequence, multichannel tunneling mechanism set a 
minimal current at a high bias, similar to a previous experimental observation. The substi- 
tution of the mispair in this case, however, will break the multichannel tunneling mechanism 
and decrease the current significantly. 
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